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Neste trabalho mostramos de forma didática a dedução da equação de onda relacionada com as chamadas 
ondas gravitacionais. Veremos que as referidas ondas surgem naturalmente como soluções das equações de Einstein 
linearizadas, após a escolha de um sistema de coordenada (gauge) adequado. Discutiremos ainda a geração de 
ondas gravitacionais e a natureza quadrupolar da radiação gravitacional emitida por um sistema binário em órbita 
circular. Por fim, aplicaremos a teoria desenvolvida nos dados observacionais da binária de Hulse-Taylor, obtendo 
assim uma noção dos valores que cercam este tipo de evento e também do quão precisa é a teoria apresentada. 
Palavras-chave: Relatividade geral, Radiação gravitacional, Binária de Hulse-Taylor. 


In this work we show a deduction of the wave equation related to the so-called gravitational waves. We will 
see that these waves arise naturally as solutions of the linearized Einstein equations, after choosing a suitable 
coordinate system (gauge). We will also discuss the generation of gravitational waves and the quadrupolar nature 
of gravitational radiation emitted by a binary system in circular orbit. Finally, we will apply the theory developed 
in the observational data of the Hulse-Taylor binary, thus obtaining notion of the values that involve this kind of 


event and how accurate the theory is. 


Keywords: General Relativity, Gravitational Radiation, Hulse-Taylor binary. 


1. Introdução 


A nossa civilização, em todas as culturas, sempre foi 
fascinada pelo conceito de espaço (céu, inferno, etc.) 
e tempo (o início, meio e fim de tudo). Ao longo do 
tempo esses conceitos se enraizaram na nossa consciência, 
dando-nos as imagens mentais que usamos. Pensamos 
no espaço como uma entidade tridimensional (largura, 
altura e espessura) que nos rodeia, e pensamos no tempo 
como uma espécie de fluxo, que guia a evolução do nosso 
universo. Juntos eles funcionam como um palco, onde se 
dá a peça que descreve o universo. Nessa peça, os atores 
são tudo que podemos imaginar - planetas, estrelas ou 
até nós mesmos. Na antiga visão Newtoniana de mundo, 
esse palco é o mesmo, seja você um ator ou um mero 
espectador. 

Porém, de acordo com Einstein em sua teoria da Re- 
latividade Geral (RG), o espaço e o tempo deixam de 
ser apenas esse palco rígido, e passam a ser um ator da 
própria peça: o chamado espaço-tempo, um contínuo de 
quatro dimensões. Esse novo “ator” é modificado pela 
presença de matéria e energia. A presença de planetas, 
estrelas e galáxias (matéria e energia) curvam a estrutura 
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do espaço-tempo de forma semelhante a uma bola pesada 
esticando uma cama elástica. 

Um dos resultados mais interessantes previsto na RG 
foi a existência de ondas gravitacionais. As ondas gravita- 
cionais podem ser pensadas, de um modo grosseiro, como 
as ondas criadas quando jogamos pedrinhas na superfície 
de um lago: elas são ondulações (perturbações) na geo- 
metria do próprio espaço-tempo geradas por qualquer 
tipo de matéria acelerada, e se propagam pelo Universo 
à velocidade da luz. Contudo, diferentemente das ondas 
no lago, as ondas gravitacionais podem se propagar no 
espaço vazio (vácuo). Além disso, uma vez que a força 
gravitacional é muito fraca, as ondas geradas são extre- 
mamente pequenas, sendo consideravelmente intensas 
apenas em regiões que concentram uma grande quanti- 
dade de massa e energia, como na região próxima de um 
buraco negro. 

O real existência dessas ondas ficou por um longo 
período restrito a discussões teóricas. Só a partir da 
década de 60 |1| os físicos passaram a buscar mais seria- 
mente a detecção dessas pequenas ondulações do espaço- 
tempo, sem ter tido nenhum resultado confiável até fe- 
vereiro do ano passado. À descoberta anunciada pelos 
detectores LIGO (Laser Interferometer for Gravitatio- 
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nal waves Observatory) Bl, e publicada na prestigiosa 
revista científica Physical Review Letters [8], não ape- 
nas confirma a existência das ondas gravitacionais como 
também fornece a primeira medição direta da existência 
de buracos negros. 

O enorme interesse gerado pela descoberta do LIGO 
trouxe a necessidade de mostrar aos alunos de graduação 
como essa descoberta se deu e também como a teoria 
de Einstein prevê essas ondas. Excelentes artigos sobre 
este tema foram publicados no fim do ano passado na 
Revista Brasileira de Ensino de Física (RBEF) e ao leitor 
sugerimos consultá-los para maiores detalhes a cerca do 
experimento e resultados afr]. 

Neste artigo nós propomos revisitar os aspectos ma- 
temáticos da RG que estão diretamente ligados ao tema. 
Nosso objetivo é apresentar de forma didática ao aluno 
de graduação o cálculo explícito das ondas gravitacio- 
nais a partir de primeiros princípios, isto é, partindo 
diretamente da equação de Einstein, acessando o regime 
linear de campo fraco e obtendo uma solução geral para 
a equação de onda resultante. Como aplicação, estuda- 
mos a radiação gravitacional emitida por um sistema 
binário de massas girando em órbita circular. Os resul- 
tados teóricos são então aplicados na análise dos dados 
do evento denominado por PSR 1913+16 ou “Binária de 
Hulse-Taylor”, cuja observação possibilitou a primeira 
detecção indireta da emissão de ondas gravitacionais [8]. 
Dessa forma, pretendemos fornecer ao leitor iniciante 
em RG uma visão concreta sobre a dedução da solução 
das equações de Einstein que dão origem às ondas gra- 
vitacionais e como usá-las para confrontar os resultados 
experimentais. Tais estudos não são geralmente cobertos 
pelos cursos de introdução a RG em nível de graduação ou 
mesmo de pós-graduação e por isso achamos conveniente 
apresentá-los aqui. 

A notação e as convenções adotadas são as mesmas da 
referência [9] e o restante do trabalho está organizada 
como segue. Na seção [2] faremos uma breve apresentação 
da equação de Einstein e e suas implicações. Na seção 
usaremos o método perturbativo para linearizar a 
equação de Einstein e mostrar que a perturbação obe- 
dece a uma equação de onda. Na seção [4] discutiremos a 
geração de ondas gravitacionais e a natureza quadrupolar 
da radiação gravitacional. Em seguida, calcularemos a 
energia emitida via radiação gravitacional por um sis- 
tema binário em órbita circular. Aplicaremos então a 
teoria desenvolvida nos dados observacionais na binária 
de Hulse e Taylor, comparando teoria e experiência. Por 
fim, na seção [5] teceremos nossas considerações finais. 


2. A equação de Einstein 


Na teoria geral da Relatividade (RG), publicada em 1916 
por Albert Einstein [10], a noção newtoniana de espaço 
como um vasto vazio no qual os eventos acontecem e, 
portanto, não tendo efeito sobre o movimento da matéria 
é transformado em uma rede “espaço-tempo” que acolhe 
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essa matéria e guia o seu curso. Nesse contexto, a própria 
ideia do que pensamos por forças gravitacionais, nada 
mais é do que uma expressão do espaço-tempo ser curvo. 
Segundo a RG de Einstein, os efeitos gravitacionais são 
incluídos em sistemas físicos através da métrica guy(x) e, 
uma vez conhecida a distribuição de matéria e energia, a 
equação que nos permite obter tal métrica é a chamada 
equação de Einstein (uma maneira interessante de se 
deduzir a equação de Einstein baseada em argumentos 
físicos de simetria pode ser encontrada na referência mj: 
T 
Guv = Tie (1) 
onde G = 6,672 x 1071N . m?/kg° é a constante gra- 
vitacional, c = 2,997 x 108m/s é a velocidade da luz 
no vácuo e Tv é o chamado tensor momento-energia. 
Trata-se do objeto mais geral que carrega a informação 
da massa, energia e momento do sistema a ser analisado. 
Do lado esquerdo da equação (1), temos o chamado ten- 
sor de Einstein (G,) que traz consigo informação acerca 
da curvatura do espaço-tempo em questão. Tal tensor é 
definido em termos de dois outros objetos 


1 
Gw = Rw = gImR, (2) 


são eles: o Tensor de Ricci 


== p= P A Pp A 
Rov = CUP = 8u — Op ToT — AURA 


(3) 
e o escalar de curvatura 
= j A A 
R= R", = g" (Ot, — 00, TRTA TAD) 
(4) 


que, por sua vez, são definidos em termos da conexão 
afim 


I$ = Lg (0,9ux T ORN z O»xJvu) . (5) 
O motivo de definirmos as expressões acima deve-se ao 
fato que posteriormente iremos linearizá-las. 

A equação de Einstein constitui um conjunto de 10 
equações diferenciais parciais não lineares em gy. Con- 
forme citado, essa equação nos permite obter a métrica 
guv quando conhecido o tensor momento-energia do sis- 
tema físico em questãåd?] Dessa forma, saberemos como o 
espaço-tempo tem sua estrutura geométrica modificada 
pela presença de matéria. Não discutiremos a inclusão 
da constante cosmológica. Considerações acerca desta 
podem ser encontradas em vários livros texto. Consulte, 
por exemplo, ou [13]. 

Vale ressaltar que as equações de campo da RG são 
consideravelmente mais complexas do que as equações 
de campo presentes na física clássica, como por exem- 
plo as equações de Maxwell do eletromagnetismo. Tal 


1Como de fato temos um conjunto de equações é comum o uso do 
termo equações de Einstein no plural. 
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complexidade surge por causa da natureza não-linear do 
campo gravitacional. Enquanto no eletromagnetismo as 
equações de Maxwell são lineares justamente pelo fato do 
campo eletromagnético não transportar carga e não inte- 
ragir consigo próprio. Em uma teoria de campo efetiva, 
o campo gravitacional é portador de energia e momento 
e, portanto, acaba gerando auto-interação, uma vez que 
no contexto da relatividade, energia e momento, além 
da massa, também são fontes de campo gravitacional. 
Uma discussão mais completa acerca desse tema pode 
ser encontrada em [9]. 


2.1. Transformação de calibre 


O tensor de Einstein, dado por p), é simétrico (pois 
Rav € guv são simétricos) e o fato de trabalharmos num 
espaço-tempo quadridimensional implica que ele possui 
apenas 10 componentes independentes, da mesma ma- 
neira que a métrica g4v. Portanto, a equação consiste 
em um conjunto de dez equações algebricamente inde- 
pendentes. Poderíamos supor que a equação de Einstein 
acrescida das condições de contorno, seria suficiente para 
que determinássemos a métrica univocamente. No en- 
tanto, isso não ocorre pois as dez componentes de G pv 
estão relacionadas pela identidade de Bianchi 


VaGt” =0, (6) 


que formam um conjunto de quatro identidades diferen- 
ciais que devem ser respeitadas. 

Como consequência, as dez equações de Einstein para 
Guv se reduzem a apenas seis verdadeiramente inde- 
pendentes, deixando-nos aparentemente impedidos de 
determinar todas as dez componentes da métrica. En- 
tretanto, se guv é solução de (1), então fazendo-se uma 
transformação de coordenada geral (x — a”) obtemos 
uma nova métrica g',, que também será solução de (1). 


Isto significa que Dis possui quatro graus de liberdade que 
dependem da escolha das coordenadas x” e não podem 
representar graus de liberdade genuinamente físicos. 

No eletromagnetismo, temos um problema semelhante 
que é o de determinar univocamente o quadripotenciaf?] 
A, através das quatro equações de Maxwell. Tal problema. 
só é resolvido quando adotamos um calibre (ou gauge?) 
particula] O problema da unicidade da métrica pode 
ser resolvida de forma análoga, basta simplesmente que 
adotemos um sistema de coordenadas específico, que 
será o nosso gauge. É a escolha desse sistema (que pode 
ser expressa por quatro condições), em adição às seis 
equações independentes que já possuímos devido a (1), 
que nos permite determinar a métrica univocamente. 
Voltaremos a esse assunto nas próximas seções, quando 
se fará necessário adotarmos um gauge. 


2 Ay = (4, A), onde ¢ é o potencial escalar e A o potencial vetor. 
3Do original em inglês. 
“Como, por exemplo, ð A! = 0 que é o chamado gauge de Lorentz. 
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3. Teoria da perturbação 


Obteremos a seguir a equação de Einstein no regime 
linear. Primeiramente, usaremos o método perturbativo 
que consiste, neste caso, em escrevermos a métrica guv 
como a de Minkowski mais uma pequena perturbação; 
substituir essa métrica perturbada na equação de Einstein 
e então linearizá-la, de modo a obter a perturbação em 
função de algum parâmetro ou classe de parâmetros. 
Dentro do regime linear, discutiremos ainda o problema 
da invariância de gauge. 


3.1. A equação de Einstein no regime linear 


A teoria linearizada da gravitação é baseada no pres- 
suposto de que em regiões do espaço-tempo distantes 
de qualquer fonte gravitacional, o campo é fraco e a 
métrica se desvia apenas levemente da de Minkowski 
(nav = (—1,1,1,1), com os demais elementos fora da 
diagonal nulos) [14]. O que torna essa suposição válida, 
é o fato de encontrarmos frequentemente na natureza 
situações nas quais a distribuição de matéria é cercada 
pelo vácuo e os corpos mais próximos estão muito afas- 
tados; de forma que o campo gravitacional pode ser tido 
como fraco pelo menos em uma região intermediária, 
onde a métrica é a minkowskiana mais uma pequena 
perturbação. Tal região é conhecida como zona de campo 
distante (ou far field) e em geral o estudo das soluções on- 
dulatórias das equações de campo são restritos às soluções 
que exibem tal região. 

Matematicamente, a aproximação mencionada acima é 
descrita como uma perturbação linear de primeira ordem 
na métrica de Minkowskf?] na forma 


Juv = Nuv + huv; (7) 
onde |h,v| << 1. A inversa da métrica fica dada por 
g” = mt” onte, (8) 


Pode-se interpretar a versão linearizada da RG (onde 
todos os efeitos são determinados em primeira ordem 
de hav), como uma teoria que descreve a propagação da 
perturbação ao longo de um espaço-tempo minkowski- 
andê] Nesse regime, as equações (5). e são escritas, 
respectivamente, como 


1 
w53 


1 
= 5 Ohya + Ouhyx = zhu), (9) 


9º (9,9ux + ON E OxJuv) 


5Esse fato também nos permitirá usar ny» e n” para levantar e 
descer índices, respectivamente. 

SEstamos considerando pequenas perturbações apenas sobre o 
espaço-tempo de Minkowski e isso não precisa ser sempre assim. Po- 
deríamos, por exemplo, considerar essas mesmas perturbações sobre 


um outro espaço-tempo qualquer de métrica g9) e, para esses casos, 


escreveríamos nossa métrica perturbada como guv = E +hpv. Ob- 
tendo assim uma teoria que descreve a propagação da perturbação 
ao longo de um espaço-tempo de métrica g9. 
necessária, por exemplo, na cosmologia. 


Essa abordagem é 
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o tensor de Riccl”] 


Rye = 9º" Ruvpo 


SM"P(0,0u ho + 9, 0,hvp — O, 0vhyap — OpOuhvo) 


1 
Ooh y + 0udoh!, — Do Oh — DO" vo), 


que como esperado é um tensor simétrico. Fazendo uma 
contração no tensor de Ricci, obtemos o escalar de cur- 
vatura 


R= g9” Rw 


1 
= z0 (08h, + 8oðpuh?, — p3h — Ohr) 


E Dan +0,0,hºt — Ə "h — Dh) 
= 3udoht” — Oh, 


(10) 


e, dessa forma, a equação de campo linearizada é escrita 
como 


30,7, + Os Bh, — Ooh — Day 
—nNuv Âp ðh + Nur Oh = 167GT pv, 


(11) 


onde Ty é o tensor momento-energia calculado na or- 
dem zero em hv, tal que sua magnitude é da mesma 
ordem da pertubação e satisfaz a lei de conservação or- 
dinária 9,7” = 0. Esta é uma suposição razoável, uma 
vez que a quantidade de energia e momento devem ser 
pequenas no limite de campos fracos [9]. Notemos que 
aqui, e no restante do trabalho, estamos assumindo c = 1 
por simplicidade. Uma vez de posse das equações de 
campo linearizadas, devemos agora lidar com a questão 
da invariância de gauge. 


3.2. Invariância de gauge 


Como explicitado anteriormente, o problema da invariância 
de gauge aparece porque a procura por uma métrica 
Juv = Nuv + huv não específica completamente o sistema 
de coordenadas no espaço-tempo; podemos ter outros 
sistemas de coordenadas em que a métrica continue a 
ser escrita como a de Minkowski acrescida de uma pe- 
quena perturbação, mas a perturbação será diferente. 
Assim, a decomposição da métrica como feita em não 
é unívoca. 

Não devemos esperar, portanto, que as equações de 
campo nos deem uma única solução, pois qualquer solução 
obtida dessas equações sempre pode gerar outras através 
de uma transformação de coordenadas (tal característica 
é assegurada pelo princípio da covariância geral). Precisa- 
mos então especificar de maneira unívoca a perturbação 
huv. Segundo [13], a transformação de coordenadas mais 
geral que preserva o limite de campo fraco é da forma 


rë — qt = x” + E(x), (12) 


TOnde h = n” hyv. 
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onde ð” /ðx” é, no máximo, da mesma ordem de mag- 
nitude que Av. 

Partindo-se da lei de transformação tensorial para a 
métrica dada por 


dx! Og!” 


“Ox? das I 


pode-se mostrar que a perturbação hyuy transforma-se, 
segundo (12), na forma 


ht = hw — 


Iuv 


e, (13) 


OM o O,» 


onde u = Nuv6”. A equação (14) é a chamada. trans- 
formação de gauge para a teoria linearizada, e a in- 
variância da teoria perante estas transformações é a cha- 
mada invariância de gauge. Pode-se provar, por inspeção 
direta, que a transformação deixa o tensor de curva- 
tura linearizado invariante. Isto é de suma importância, 
uma vez que a física do problema está principalmente 
no tensor de curvatura (ou de Riemann) R,»po: Conse- 
quentemente, a equação de campo dada por (11) também 
é invariante sobre a mesma transformação” Assim, fi- 
xando as quatro funções representadas por £”, estaremos 
especificando univocamente a perturbação. 


(14) 


4. Radiação gravitacional 


Nesta seção discutiremos a geração de ondas gravitacio- 
nais e a natureza quadrupolar da radiação gravitacional. 
Por fim, aplicaremos a teoria aqui desenvolvida aos dados 
observacionais da binária de Hulse e Taylor. 


4.1. Produção de ondas gravitacionais 


Conforme mencionamos acima, a equação de Einstein 
linearizada é invariante de gauge e sua forma parece 
proibitiva para qualquer tentativa de solucioná-la. En- 
tretanto, a liberdade de gauge nos permite escolher um 
sistema de coordenada que possa simplificar bastante a 
forma da equação. Uma das escolhas mais comuns é o 
chamado sistema de coordenada harmônico (ou simples- 
mente gauge harmônico), definido por 


id =0, (15) 
cuja forma linearizada é dada por 
i 1 
o huv = gov. (16) 


Tal escolha particular de gauge é sempre possível pois 
mesmo que a perturbação original h,, não a obedeça, 


SEssa invariância a qual nos referimos é semelhante a invariância de 
gauge do eletromagnetismo. Na teoria eletromagnética, a física está 
contida nos campos eletromagnéticos, que são geralmente escritos 
em termos do quadripotencial A”. O problema reside no fato de 
o quadripotencial não ser unívoco (diferentes quadripotenciais 
podem gerar os mesmos campos). Tal problema é resolvido fixando 


o quadripotencial através de uma transformação, geralmente A — 
Ap +ôpà. 
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podemos usar (14) e obter uma nova perturbação huv 
que satisfaz a condição (16). De fato, basta realizar uma 
transformação de coordenada que seja solução de 


(17) 


1 
Ev = 9“ hyv — garh- 


Assim, no gauge de harmônico, a equação de Einstein 
(11) assume a forma 


(18) 


1 
2 no Zh) -= —167GT pv. 


A equação anterior indica que seja conveniente intro- 
duzirmos a perturbação de traço reversd?] definida por 


= 1 
huv == huv = gevh (19) 
sendo esse nome é justificado pelo fato de que h = —h. 


Para tal perturbação, a condição de gauge harmônico 
implica em “hy = 0. 

Devido à enorme semelhança com o gauge de Lorentz 
(ð A” = 0) do eletromagnetismo, a condição acima 
(ath, = 0) recebe o mesmo nomd! Assim, no gauge 
de Lorentz, a equação de Einstein assume a forma sim- 
plificada 


av = —161GT pv, (20) 


sendo portanto uma equação de onda para cada uma das 
componentes da perturbação. 

Dada uma fonte T,,, nosso objetivo é encontrar a 
solução geral para a equação e trabalhar as proprie- 
dades da radiação emitida por um sistema em diferentes 
tipos de movimento. Matematicamente, isso é semelhante 
ao problema de encontrar a solução para as equações de 
Maxwell na presença de fontes, onde tínhamos de resol- 
ver ODA” = —4r J” para determinar o quadripotencial. A 
solução da equação é obtida a partir do método da 
função de Green, cuja a solução é 


= 1 
Faltz) =4G f Tult- -yl y)dy, (21) 


onde t = x°. A quantidade t, = t — |x — y| é chamada de 
tempo retardado e sua ocorrência no argumento de 
implica que a perturbação gravitacional se propague com 
velocidade unitária, que assumimos ser precisamente a 
velocidade da luz no vácuo. 

A solução pode ser interpretada da seguinte ma- 
neira: a perturbação no campo gravitacional em (t,x) é 
a soma das influências da energia e momento da fonte 
no ponto (tr,x — y) no cone de luz passado. Percebe- 
mos que enquanto a gravitação newtoniana possui uma 
dinâmica trivial (isto é, todas as interações acorrem de 


9Podemos reconstruir a perturbação original a partir de Pis Assim, 
nenhuma informação é perdida ao fazermos essa mudança. 
10 Vale ressaltar que a perturbação original não é transversa nesse 
gauge pois Oh!” = 50h. 
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forma instantânea), a gravitação einsteiniana é não tri- 
vial, ou seja, a perturbação gerada na fonte deve levar um 
certo tempo até chegar a outro ponto do espaço-tempo. 
Vale ressaltar que existe ainda uma outra solução para a 
equação que está relacionada com o tempo avançado 
ta = t+ |x — y|, no entanto a descartamos, pois viola 
causalidade. 


4.2. A natureza quadrupolar da radiação 
gravitacional 


Nosso próximo passo será extrair informações acerca 
da radiação gravitacional emitida por uma fonte Tyv. 
Antes de prosseguir, será conveniente eliminarmos a de- 
pendência de Ty no tempo retardado, o que pode ser 
conseguido através de uma transformada de Fourier na 
variável tempo [12]. Se Ñ, (w, y) e huv(w,x) são, res- 
pectivamente, as transformadas de Fourier de Tyu (t, y) e 


hyv(t,x), então 


1 ds 
huv(tw,x) = ar far Ru (É, Xx) 


= 4G dt ot fayt be — y|,Y) 
v2m |x -y| 
4G | civla- 

a di ea ml 


elulx—yl] 
= 16 [ay Tulio, y). (22) 
bx -yl 
Na equação acima, a primeira linha representa simples- 
mente a definição da transformada de Fourier no tempo, 
na segunda linha nós inserimos a solução (21), na ter- 
ceira linha usamos a definição do tempo retardado t, = 
t—|x— y| e o fato de dt, = dt para implementar uma 
mudança na variável de integração de t para t,, e final- 
mente, na quarta linha usamos novamente a definição da 
transformada de Fourier. 

Vamos agora nos restringir a situação na qual a fonte 
encontra-se isolada, muito distante e é composta apenas 
por matéria não relativista] Tais aproximações são 
adequadas e suficientes para o estudo da radiação emitida 
por um sistema binário em orbita circular que trataremos 
a seguir. 

Pondo a origem do sistema de coordenadas perto da 
fonte, a aproximação adotada nos permite trocar na 
equação |x — y| por r = |x|, o que faz com que o 
termo ei *=Y| /|x — y| possa ser reescrito como ei” /r. A 
aproximação de baixas velocidades (quando comparadas 
à velocidade da luz) nos permite assumir que r é prati- 
camente constante, pois demora muito tempo até que a 
distância entre a fonte e o observador tenha alguma dife- 
rença significativa. Com tudo isso, podemos reescrever 
(22) como 


iwr 


hur(w,x) = 16º [ay Turvo, y). (23) 
r 


11 Matéria cuja a velocidade típica é muito menor que a da luz. 
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Felizmente, precisaremos calcular apenas as componen- 
tes espaciais de h,y(w,xX), e as razões para isso são duas: 
uma de caráter físico e outra matemático. A motivação 
de caráter físico reside no fato de que a lei de conservação 
O, Tt” = 0 implica que todas as integrais sobre Tº” são 
conservadas e independentes do tempd?] Dito de outra 
maneira, a transformada de Fourier de T°” com respeito 
à variável tempo conterá apenas termos estáticos (w = 0). 
Isso nos mostra que a parte dependente do tempo da 
métrica, que é o que caracteriza a radiação, está contida 
nas componentes espaciais da perturbação [12]. 

A razão de ordem matemática está associada com 
o gauge adotado ð h” = 0 (que implica 9,7” = 0 
por causa de (20)). Esta por sua vez, estabelece uma 
relação entre as componentes da perturbação no espaço 
de Fourier, dada por: 


iu (00,%) + O (00,3) = 0. (24) 


03 
Tal equação nos mostra queh pode ser determinado a 


partir de nº . Além disso, juntando o conhecimento de que 
di = A com a equação (24), é possível determinar no 
Portanto, precisamos calcular apenas as componentes 
espaciais da perturbação, pois todas as outras podem ser 
obtidas a partir delas. 

Dadas as devidas motivações, passamos então a resolver 
(23) apenas para as componentes espaciais. Começamos 
por efetuar uma integração por partes 


fè T” (w,y) = fë Ər (y TS) 
— fë y’ (aL TI) a (25) 


onde o primeiro termo do lado direito da igualdade é uma 
integral de superfície que irá desaparecer pelo fato de 
a fonte se encontrar isolada; enquanto o segundo termo 
pode ser relacionado com Ti, uma vez que a lei de 
conservação ð T!” = 0, quando escrita no espaço de 
Fourier, impõe que 


Op Tt = iw T. (26) 
Assim, 
fë T“ lw,y) = -iw [ay y To 
iw da as 
es dê ipoj j TOi 
a / y (yT ei) 
dy dê [ou (y'yi TO) — ini (O, Ok 
mi y [de(y'y yy (Op TO] 
w? onga 
a fë yiyi T. (27) 
12 Bxplicitamente: | ato = f d8zooTo” = — f ProT” = 


— $ ds;T” = 0, pois as fontes são localizadas. 
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A segunda linha é justificada por conta da simetria dos 
índices à e j, enquanto as duas últimas foram novamente 
aplicações da integração por partes e a equação (26). É 
convencional na literatura [9/12/13], definir o chamado 
tensor momento de quadrupolo da distribuição de matéria, 
como sendo 


I (t) = fë gy Tt, y), (28) 


que juntamente com (27) nos permite escrever a equação 
(23) na forma 


iwr 


= e = 
hij(w,x) = —2Gw? Ii; lw). (29) 

Retornando para a dependência temporal, a equação 
anterior é reescrita como 


z= 2G d? 
hij(t, x) = + gg) (30) 


que é a chamada fórmula de quadrupolo onde tr = t — r. 

Esse resultado nos revela que a radiação gravitacional 
emitida por uma fonte, nas condições descritas, é propor- 
cional à derivada temporal de segunda ordem do tensor 
momento de quadrupolo da distribuição de matéria. Po- 
demos observar que isso contrasta com a radiação eletro- 
magnética, em que a maior contribuição é proporcional 
à derivada temporal de segunda ordem do momento de 
dipolo da distribuição de cargas, e o motivo disso é fácil 
de entendermos. 

O momento de dipolo no eletromagnetismo é definido 
usando-se a distribuição de cargas, e para um sistema 
feito de cargas com diferentes taxas q/m, ele pode mudar. 
No caso gravitacional, o momento de dipolo d de um 
sistema de massas é dado por: 


d = MRom, (31) 


onde M é a massa total do sistema e Rcm é o vetor que 
localiza o centro de massa. No caso de distribuições isola- 
das, temos que o centro de massa se move com velocidade 
constante e, portanto, Rom = 0. Concluímos, então, que 
não podemos ter radiação gravitacional proveniente do 
momento de dipolo [12]. Assim, na decomposição mul- 
tipolar, o termo de ordem mais baixa para a radiação 
gravitacional é o momento de quadrupolo. Vale ressaltar 
também que o momento de quadrupolo é geralmente 
menor que o de dipolo, e é por essa razão, bem como o 
fraco acoplamento de matéria a gravidade, que implica 
que a radiação gravitacional é tipicamente muito menor 
do que a eletromagnética [9]. 


4.3. A energia transportada por ondas 
gravitacionais 


Vamos agora a uma das questões mais delicadas da te- 
oria: a energia transportada por ondas gravitacionais. 
Assunto por sinal muito importante, pois nas teorias de 
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campo tratadas por intermédio da RG, temos que o fator 
decretório para a existência de ondas não é apenas que 
os campos sejam dependentes do tempo; mas que, além 
disso, haja transporte de energia e momento sem haver 
transporte de matéria. No entanto, essa é uma questão 
muito difícil de se lidar [L5|[16]. 

As complexidades desse assunto não são apenas de 
origem técnica, mas também conceituais, uma vez que 
o próprio princípio da equivalência - pilar básico da 
RG - não nos permite definir a energia associada a um 
campo gravitaciona[B] Neste trabalho, não discutiremos 
o mecanismo pelo qual calcula-se a energia transportada 
por ondas gravitacionais. Apenas usaremos a relação, já 
existente na literatura [9/12], de que a potência irradiada 
(Lew) por sistemas físicos devido a emissão de ondas 
gravitacionais é dada por 


_dEew  G/EIS SIA, 
vu gia qt A 


onde TZ” é a versão transversal e de traço nulo do mo- 
mento de quadrupolo calculado, como antes, no tempo 
retardado. A fórmula acima representa a média temporal 
da taxa de energia que é trocada, e fontes irradiantes 
estão perdendo energia, por isso o sinal negativo na igual- 
dade. 


4.4. Ondas gravitacionais emitidas por um 
sistema binário em órbita circular 


Como aplicação dos resultados obtidos até agora, usare- 
mos a fórmula para calcular a energia emitida, via 
radiação gravitacional, por um sistema bem simples: dois 
corpos de massas mı e ma efetuando órbitas circulares 
em torno de um centro de massa em comum. Tal situação 
é descrita na figura [I] 

Faremos uma análise dessa situação através da mecânica 
newtoniana, válida em baixas ordens. Neste caso, os 
parâmetros orbitais são completamente determinados 
pela massa reduzida u e pela distância ly entre os dois 
corpos envolvidos. Baseando-se na figura acima para a 
nossa análise, onde adotamos a origem do sistema de 
coordenadas coincidindo com o centro de massa do sis- 
tema, podemos escrever, em termos da massa reduzida 
u = mima/M e a separação orbital lo = rı + r2, que as 
distâncias rı e r2 valem respectivamente 


mal, 

r= E (33) 
mal 

r9= TE (34) 


onde M = mı + ms é a massa total do sistema. 


13Na RG a definição de energia (componente 00 do tensor momento- 
energia) depende do observador. Uma forma simples de ver isto 
é a seguinte: não sabemos dizer se a gravidade sentida é devido 
ao referencial está acelerado ou por existir realmente um campo 
gravitacional. 
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Figura 1: Dois corpos de massas mı e mz efetuando órbitas 
circulares em torno de um centro de massa em comum. 


A frequência orbital Q é obtida igualando-se a força 
gravitacional newtoniana, a resultante centrípeta. Para 
cada uma das massas teremos 


Mama 2 Molo 
Ea O) 
G R ma M , (35) 
Mama 2 malo 
G—— = MN 36 
Iê ma M , ( ) 
o que implica que para ambas as massas 
GM 
0 
e, portanto, as equações de movimento serão 
mal mal 
zı = a cos(Qt) , y= o sen(Mt), (38) 
l l 
To= r cos(Qt) , yYy2= TA sen(Qt). (39) 


De posse dessas equações estamos aptos a computar o 
tensor momento de quadrupolo, definido por (28). Para 
este caso temos 


2 
TO = SD mudla — co)ôly — yn)ô(2), (40) 
n=1 
o que resulta em 
Las mi [a2ô(a — xı)ô(y — yr )ó(z)dadydz 


+ mo I v2ô(a — z2)ô(y — y2)ô(z)drdydz 
= mir? + mz? 

= pulĝcos?’ (Qt) 

- 5 cos(20t) + 56 (41) 
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Aplicando-se o mesmo procedimento para as outras com- 
ponentes, encontramos que todas são nulas, exceto 


Lu —Slgcos(204) = 5i (42) 


Tey = Slgsen(20). (43) 


No entanto, note que precisaremos apenas da parte 
temporal dessas componentes e, portanto, as contas serão 
efetuadas tomando-se 


Tss = —Iyy = Siôcos(201), (44) 


Ly = 5 sen(29t). (45) 


Antes de prosseguirmos é interessante notar que a 
frequência das ondas gravitacionais emitidas é o dobro 
da frequência orbital. Em uma forma matricial o tensor 
momento de quadrupolo é escrito como 


cos(2Qt) sen(2Qt) O 
lj= 5 sen(20t) —cos(2Qt) 0 |, 
0 0 0 


ou de forma mais reduzida 


H 
Ii; = 510415» (46) 
onde A;; são os elementos da matriz apresentada acima. 
Substituindo esse resultado em (32), encontramos que a 
magnitude da radiação gravitacional emitida pelo sistema 
binário é dada por 
L _ dEgw _ 32 Gu M? 
Word 5 el 

onde reinserimos a constante c no resultado acima. 

A ordem de magnitude da amplitude dessas ondas 
pode ser calculada substituindo-se (46) na fórmula do 
quadrupolo, o que resulta em 


_ 4uMG? 
B rloct i 


(47) 


ho (48) 


4.5. A binária de Hulse e Taylor 


Em 1974, no Observatório de Arecibo, dois astrônomos 
Russell Hulse e Joseph Taylor fizeram uma importante 
descoberta: uma estrela binária formada por um pulsar 
e outra estrela de nêutrons que orbitam um baricentro 
comum. Tal sistema foi designado por PSR 1913 + 16, 
e foi o primeiro a possibilitar uma detecção indireta de 
ondas gravitacionais. O que faremos agora é aplicar os 
resultados obtidos na subseção anterior a este sistema, 
a fim de termos uma ideia dos valores envolvidos nesse 
tipo de evento; por exemplo, a taxa da liberação de 
energia devido à emissão de ondas gravitacionais. Como 
consultado em [8/17], os dados do referido evento são: 


my  1,441M5, 
ma ~ 1,387M6, 
lo = 0,19 - 101° m, 
P = Th 45m Ts, 
r= 1,5-10% m, 
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onde Mo = 1,99 x 10%kg e Ro = 6,96 x 108m são, 
respectivamente, a massa e o raio do Sol. Os outros 
parâmetros acima são os mesmos descritos na subseção 
anterior; P o período e r a distância binária-Terra. De 
posse desses dados, podemos concluir que a frequência 
vew das ondas emitidas nesse evento são da ordem de 


29) 


vew = —  7,2-:10"º Hz, (49) 
27 


além disso, usando-se (47) e (48) encontramos também 


Lew = 7,0-10W, 
ho ~ 6,12: 10728. 


Os resultados acima nos permitem tirar algumas con- 
clusões imediatas. Primeiro, a amplitude de tais ondas é 
ínfima. Incrivelmente, por mais alta que seja a energia 
irradiada nesse evento, nem mesmo os observatórios do 
Advanced LIGO seriam capazes de detectar essas ondas. 
As curvas de sensibilidade do LISA e LIGO são mostradas 
na Figura 2] Neste gráfico, a amplitude adimensional ca- 
racterística está em unidades “por raiz quadrada de Hz”. 
Aplicando os valores obtidos em e na fórmula 
de conversão ho vzs, obtemos o valor aproximado de 
7,2 x 10-21Hz7 2 para o sistema binário PSR 1913 + 16 
( para maiores detalhes sobre as convenções utilizadas na 
literatura para as curvas de sensibilidade de ondas gravi- 
tacionais, sugerimos ao leitor consultar a referência [18)). 
Notemos que sua localização está fora das regiões de 
detecção dos aparelhos [4] Entretanto, outros sistemas 
binários mais próximos da Terra poderiam ser detecta- 
dos. Por exemplo, os sistemas detectados pelo LIGO, 
teriam sido detectados pelo LISA, uns 10 anos antes, se 
ele já estivesse existindo (ver seminário apresentado na 
conferência Amaldil2 [19], recentemente). 

Observamos ainda que a perda de energia devido à 
emissão de ondas gravitacionais é altíssima. Para se ter 
uma ideia de quanta energia estamos falando, o próprio 
Sol irradia uma quantidade equivalente a 4,3: 10% W 
através de radiação eletromagnética, mas apenas 5:103 W 
via radiação gravitacional. É de se esperar que com a 
perda de toda essa energia as componentes da binária 
acabem se aproximando, o que acarreta na diminuição 
do período orbital. Para finalizar o trabalho, o que fare- 
mos agora é aproveitar todos os resultados obtidos até 
aqui para calcular essa diminuição no período. Isso nos 
possibilitará ter uma noção de quando irá acontecer a 
coalescência das referidas componentes. 

A expressão foi obtida a partir de (82). Isto implica 
que ela deve ser considerada como uma média de muitos 
comprimentos de onda ou, equivalentemente, sobre um 
número muito grande de períodos. Portanto, a fim de 
que Lew seja uma grandeza definida, devemos estar em 


Podemos criar nossa própria curva de sensibilidade e as regiões 
de localização dos eventos através dos aplicativos online dis- 


poníveis em http://www.srl.caltech.edu/-shane/sensitivity/ 
MakeCurve.html e http://rhcole.com/apps/GWplotter/. 
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Figura 2: Em vermelho, temos a curva de sensibilidade dos equi- 
pamentos LISA (Laser Interferometer Space Antenna Project). 
Em azul, a curva de sensibilidade dos equipamentos do Advanced 
LIGO. Apenas eventos localizados acima de tais curvas podem 
ser detectados. No entanto, como pode ser visto, PSR 1913+16 
não se encontra em nenhuma dessas regiões e, portanto, não 
poderia ser detectado. 


uma região onde os parâmetros orbitais não mudam sig- 
nificativamente durante o intervalo de tempo necessário 
para executar tal média. Esta hipótese é chamada de 
aproximação adiabática e, certamente, é aplicável a siste- 
mas tais como PSR 1913 + 16 que estão muito longe da 
coalescência5] 

No regime adiabático, o sistema tem tempo de ajustar 
a órbita para compensar a energia perdida pela emissão 
de ondas gravitacionais com uma mudança na energia 
orbital, de tal maneira que: 


dEorb 
dt 


Vejamos agora as consequências dessa equação. 

A energia orbital é igual à energia cinética Æe de cada 
uma das componentes, acrescida da energia potencial U 
de ligação (que neste caso é de origem gravitacional), isto 
é 


ciano; (52) 


Eorb = Ec Sr U, (53) 
onde 
1 2 1 2 
Ec = au (Qr1) + 22 (Qro) 
ola mamalg mâmalo 
2 M2 M? 
GuM 
4 
= (54) 
E G GuM 
HE mama Z H (55) 
lo lo 
Substituindo e em (53), obtemos 
GuM 
E pila (56) 


2lo * 


150s dados mostrados acima indicam que essa coalescência está 
longe de ocorrer. O que queremos é ter ideia do quão longe. 
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e tomando-se a derivada temporal da expressão acima 


dEor GuM (1 dl 1d! 

5 = E x = Eorb z $ (57) 
dt 2lo lo dt lo dt 
No entanto, a expressão acima pode ser escrita de uma 


maneira mais sugestiva, para isso, reescreveremos dlo/dt 
em termos da derivada temporal de Q. De (37), temos 


1 dQ 3 dlo 
Qd 2 dt Es) 
e substituindo em (57), ficamos com 
dEor 2 Eorp dQ 
bo b (59) 


di 3 Q dt’ 


Por outro lado, a frequência orbital Q é relacionada com 
o período P pela relação Q = 27r P, o que implica em 


1 dQ 1 dP 
= 60 
Q dt P dt’ o 
pondo isso na equação (59), obtemos 
dP 3 P dE 
= 1 
dt 2 Eowp dt (61) 
ou, fazendo uso de 
dP 3 P 
— = -—— Lev. 2 
dt ama C” (62) 


A relação acima nos permite calcular o quanto o 
período da órbita muda devido a emissão de ondas gra- 
vitacionais. Supondo que nossa órbita seja circular e, 
portanto, que os dados acima são verdadeiros, temos que 


P = 27907 s 
Ex» ~ —1,4 -10% J 
Lew =7,0 -10 W 
logo 


dP 
—— ~ —2, 2. 107"s/s. 
q : 0" ““s/s (63) 


No entanto, o valor de dP/dt foi obtido a partir de 
medições. Tal valor, medido com uma precisão altíssima, 
foi encontrado depois de aproximadamente três décadas 
de monitoramento e vale: 


P 
= = — (2,4184 + 0, 0009) 10-12 s/s, (64) 


o que implica em uma taxa de redução de período orbital 
de aproximadamente 76,5 microssegundos por ano. Com 
isso, o tempo de vida calculado para a espiral final é de 
300 milhões de anos [8]. 

Como pode ser visto, existe uma grande diferença entre 
o valor obtido teoricamente e o medido. Na verdade, tal 
discrepância já era de ser esperada. Isto porque a órbita 
real do sistema não é circular, mas sim uma elipse de 
excentricidade e ~ 0,617 [8]. Como explicado em [12], se 
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refizéssemos os cálculos usando as equações de movimento 
apropriadas para uma órbita excêntrica, obteríamos 


dP 
nv —2,4. 1071s 
T i 0" ““s/s (65) 


que é incrivelmente próximo do valor medido. 

Em 2004, Joseph H. Taylor e Joel M. Weisberg publi- 
caram uma nova análise dos dados experimentais, con- 
cluindo que a disparidade de 0,2% entre os dados e 
os resultados previstos é devido a constantes galácticas 
pouco conhecidas, incluindo a distância do Sol ao centro 
da galáxia, movimento próprio do pulsar e sua distância 


da Terra r 


5. Conclusão 


O objetivo deste trabalho foi fazer um estudo introdutório 
sobre ondas gravitacionais. Começamos por fazer, na 
seção [2] uma breve revisão sobre a equação de Einstein 
da RG. Em seguida, na seção |3| aplicamos o método 
perturbativo nas equações de Einstein e em seguida as 
linearizamos. Na seção |4| estudamos a geração de ondas 
gravitacionais. Para isso, adotamos o gauge harmônico 
o que nos permitiu reescrever a equação de Einstein na 
forma de uma equação de onda com termo de fonte, isto 
é, considerando a presença de matéria. Passamos então 
a discutir a natureza quadrupolar da radiação e suas di- 
ferenças em relação à radiação eletromagnética que é de 
natureza dipolar. Por fim, computamos a energia emitida 
via radiação gravitacional por sistemas binários em órbita 
circular e aplicamos nossos resultados aos dados experi- 
mentais da PSR 1913 + 16 (Binária de Hulse e Taylor). 
Os dados desse pulsar binário levam à previsão teórica de 
que a potência irradiada nesse evento, devido exclusiva- 
mente às ondas gravitacionais, é da ordem de 7: 102) W. 
Para se ter uma ideia do que significa essa quantidade de 
energia, o próprio Sol irradia uma quantidade equivalente 
a 4,3-1026 W através de radiação eletromagnética, mas 
apenas 5 - 103 W via radiação gravitacional. 

Outro dado de suma importância obtido aqui foi a 
diminuição do período orbital dessa binária, algo em 
torno de dP/dt = —2,2-10"13s/s. No entanto, este re- 
sultado não é tão próximo do valor medido dP/dt = 
—(2,4184 + 0, 0009) - 10712s/s (uma redução de aproxi- 
madamente 76,5 us por ano). Tal medição foi resultado 
de aproximadamente três décadas de muitas observações 
e realizada com uma precisão altíssima. De fato, esse 
erro já era esperado, pois a órbita do sistema real não 
é circular e sim elíptica, com excentricidade e œ 0,617. 
Como pode ser visto em [12], fazendo-se novamente os 
cálculos, mas dessa vez usando as equações de movimento 
apropriadas para uma órbita excêntrica, iríamos encon- 
trar dP/dt = —2,4-10-12s/s, o que concorda com o 
resultado experimental. 
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